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ABSTRACT 

We report the discovery of an extremely broad Ha emission line in the LINER nucleus of NGC 4579. 
From ground-based observations, the galaxy was previously known to contain a Type 1 nucleus with a 
broad Ha hue of FWHM « 2300 km s'^ and FWZI « 5000 km s'^ New spectra obtained with the 
Hubble Space Telescope (HST) and a O'.'2-wide slit reveal an Ha component with a FWZI of ~ 18,000 
km s~^. The line is not obviously double-peaked, but it does possess "shoulders" on the red and blue 
sides which resemble the Ha profiles of double-peaked emitters such as NGC 4203 and NGC 4450. This 
similarity suggests that the very broad Ha profile in NGC 4579 may represent emission from an accretion 
disk. Three such objects have been found recently in two HST programs which have targeted a total 
of 30 galaxies, demonstrating that double-peaked or extremely broad-line emission in LINERs must be 
much more common than would be inferred from ground-based surveys. The ratio of the narrow [S II] 
AA6716, 6731 lines shows a pronounced gradient indicating a steep rise in density toward the nucleus. 
The direct detection of a density gradient within the inner arcsecond of the narrow-line region confirms 
expectations from previous observations of lincwidth-critical density correlations in several LINERs. 

Subject headings: galaxies: active — galaxies: individual (NGC 4579) — galaxies: nuclei — galaxies: 
Seyfert 



1. INTRODUCTION 

In the twenty years since LINERs were first identified as 
a major component of the extragalactic population (Heck- 
man 1980), most research on these objects has been di- 
rected toward resolving the question of whether they rep- 
resent a type of accretion-powered active galactic nucleus 
(AGN), or whether they are instead powered by purely 
stellar phenomena. The LINER class may be heteroge- 
neous, but it is now very clear that at least some LIN- 
ERs must be genuine AGNs (for a review, see Ho 1999a). 
An exciting recent development has been the detection of 
very broad, double-peaked Ha emission lines in a handful 
of nearby LINERs (Storchi-Bergmann, Baldwin, & Wilson 
1993; Bower et al. 1996; Shields et al. 2000; Ho et al. 2000), 
a phenomenon which is most readily interpreted in terms 
of emission from a relativistic accretion disk. 

NGC 4579 is a Virgo cluster Sb-type spiral with a red- 
shift of cz = 1520 km s~^ and a distance of 21 Mpc ac- 
cording to recent TuUy-Fisher measurements (Gavazzi et 
al. 1999). It has long been known to host a LINER nucleus 
with broad wings on the Ha emission-line profile (Stauf- 
fer 1982; Keel 1983; Filippenko & Sargent 1985). Using 
a ground-based spectrum obtained at the Palomar 5-m 
telescope. Ho et al. (1997b) found that the broad compo- 
nent of Ha has a full-width at half maximum (FWHM) 
of 2300 km s'^. Hubble Space Telescope (HST) Faint Ob- 
ject Spectrograph data have revealed substantially broader 
permitted lines in the ultraviolet (UV): for example, C IV 



A1549 has FWHM = 6600 km s'^ and full-width at zero 
intensity (FWZI) of 13,000 km s^^ (Barth et al. 1996). In 
this paper, we describe an HST optical spectrum of the 
NGC 4579 nucleus, which provides new insights into the 
properties of this low-luminosity AGN. Our analysis of the 
rather complex emission-line kinematics will be presented 
elsewhere. 

2. OBSERVATIONS 

We observed NGC 4579 with the Space Telescope Imag- 
ing Spectrograph (STIS) on 1999 April 21 UT as part of a 
program to study the emission-line kinematics of a sample 
of six nearby AGNs. The galaxy was observed at five par- 
allel positions with the STIS O'.'2-wide slit, with a gap of 
O'.'OS between adjacent slit positions. Six spacecraft orbits 
were divided roughly equally among the five slit positions, 
and the exposure time for the central slit position was 
2700 s. The telescope was oriented with the STIS slit at a 
position angle of 95?3. 

The observations were obtained with the G750M grat- 
ing and a tilt setting of 6300-6860 A, giving coverage of 
[O I] AA6300, 6363, the Ha-h[N II] emission blend, and the 
[S II] AA6716, 6731 doublet. For this setting, the point- 
source spectral resolution is 0.9 A. The CCD was read 
out in unbinned mode, for a pixel scale of 0.554 A along 
the dispersion direction and 0'.'0507 along the spatial axis. 
Three separate exposures were combined with rejection of 
cosmic-ray hits. The data were calibrated by the stan- 
dard pipeline at STScI. To remove remaining hot pixels. 
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we performed one additional step of cleaning the data us- 
ing a filtering algorithm which replaced deviant pixels with 
a local median value. The cleaned images were then geo- 
metrically rectified using the STIS calibration software. 

3. BROAD Ha EMISSION 
3.1. Observed Properties 

Figure 1 shows the HST spectrum of the NGC 4579 nu- 
cleus. The extraction is centered on the nucleus and is 
a sum of four CCD rows (0'.'2 total) in the 2-dimensional 
spectrum. The background level (including starlight from 
the galaxy bulge) is negligible in comparison with the nu- 
clear brightness, and no background subtraction has been 
applied. In addition to the usual narrow-line features ex- 
pected in a LINER spectrum, there is a surprisingly broad 
and asymmetric component to the Ha emission. On the 
blue side, the profile has a "shoulder" at 6520 A, just blue- 
ward of the base of the [N II] A6548 line. The red side of 
the profile does not exhibit such well-defined features, al- 
though there is a shallow shoulder on the profile at ~ 6670 
A. These shoulders are real features and not merely arti- 
facts of the data reduction or geometric rectification, as 
they also appear in spectral extractions of the raw data, 
even for large extraction widths. The blue edge of the line 
appears to be at ~ 6420 A rest wavelength, while its red 
tail extends to the end of the spectrum at 6820 A, giving 
FWZI « 18,000 km s~^. For comparison. Figure 1 also 
displays the ground-based spectrum from Ho, Filippenko, 
& Sargent (1997a), which was obtained at the Palomar 
5-m telescope on 1984 February 12 UT. 

We performed profile fits to remove the narrow-line con- 
tributions of Ha, [N II], and [S II]. To fit the narrow lines, 
we first fit a spline to the approximate shape of the very 
broad Ha emission profile, and subtracted this component 
from the spectrum, leaving only the narrow lines and a 
possible contribution from a "normal" broad Ha compo- 
nent. The [S II] doublet was fitted with a pair of Gaus- 
sians. For the Ha-|-[N II] blend, we used a four-Gaussian 
fit representing narrow and broad Ha and the two [N II] 
lines, which were constrained to have the same velocity 
width and a height ratio of 3:1. The narrow Ha, [N II], 
and [S II] hues have FWHM = 520-570 km s-\ while the 
normal broad component of Ha has FWHM = 2100 km 
s~^, very similar to the value of 2300 km s~^ found by Ho 
et al. (1997a). 

After subtracting the fitted components, the remaining 
very broad component of Ha has / — 1.1 x 10~^^ erg cm~^ 
s~^, corresponding to a luminosity of 5.8 x 10^^ erg s~^ 
for a distance of 21 Mpc. Since there have been no pre- 
vious optical HST spectra of the NGC 4579 nucleus, it is 
not possible to determine whether this component has al- 
ways been present or whether it is a transient feature. The 
amplitude of the broad feature is sufficiently low that it 
could be present but undetected in the Palomar spectrum. 
In general, it is extremely difficult to detect such broad, 
low-amplitude features in ground-based spectra, because 
of the overwhelming contamination by starlight from the 
bulge of the host galaxy. To make matters worse, ellip- 
tical galaxies and spiral bulges generally display a broad 
continuum hump in the region 6400-6700 A (e.g., Filip- 
penko & Sargent 1985). The starlight subtraction method 
employed by Ho et al. (1997a) was designed to yield zero 



continuum flux in the subtracted spectrum, and this pro- 
cedure could have removed a faint, very broad Ha feature 
by oversubtracting the stellar continuum template. 

No other broad optical lines have been detected in NGC 
4579, but the HST UV spectrum (taken in 1994 Decem- 
ber) clearly demonstrates that other permitted lines are 
very broad as well. The UV spectrum has low signal-to- 
noise ratio (S/N) in the grating setting that covers Ly a 
and C IV, but emission from background starlight in this 
wavelength range is negligible. As shown in Figure 2, the 
Ly a and C IV profiles have blue wings which extend to 
~ 6000 km s~^ from the systemic velocity, comparable to 
the maximum extent of the blue wing of Ha. The red wing 
of Ha extends to velocities roughly 4000 km s~^ greater 
than the red wing of C IV. It is difficult to discern the 
maximum extent of the red wing of Ly a due to the low 
S/N and a possible contribution from N V A1240, but Ly 
a emission is seen out to perhaps 9000 km s~^ from the 
systemic velocity. 

3.2. Relationship to the Central Engine 

Double-peaked broad Ha emission lines have previously 
been detected in six LINERs as well as in several radio 
galaxies and radio- loud quasars; see Eracleous & Halpern 
(1994), Filippenko (1996), and Ho et al. (2000) for reviews 
of their properties. Models of emission from rotating, Ke- 
plerian accretion disks provide the best explanation for the 
origin of the double-peaked Ha lines in these objects (e.g., 
Halpern & Filippenko 1988; Chen, Halpern, & Filippenko 
1989; Halpern et al. 1996; Storchi-Bergmann et al. 1997). 
The width of the broad-line profile in NGC 4579 is com- 
parable to the linewidths in these double-peaked emitters, 
which have typical FWZI « 15,000-20,000 km s'^. While 
the NGC 4579 Ha profile is not clearly double-peaked, the 
shoulders at —2500 and 5000 km s~^ are very similar to 
the features detected in NGC 4203 (Shields et al. 2000) 
and NGC 4450 (Ho et al. 2000). The overall similarities in 
the line profiles suggest a common origin for the broad-line 
emission in these three galaxies. Their profile shoulders are 
much less pronounced than the broad humps in "classic" 
double-peaked emitters such as Arp 102b (e.g., Halpern & 
Filippenko 1988), but Ho et al. (2000) have shown that a 
relativistic accretion disk model can fit the double-peaked 
profile in NGC 4450. 

Assuming that the very broad emission from NGC 4579 
originates in an accretion disk, the inner radius of the line- 
emitting portion of the disk can be estimated from the 
maximum velocity extent of the line. For a disk inclination 
of 45°, the projected velocity of the red wing of Ha (12,000 
km s^^) corresponds to a radius of 160 Schwarzschild radii 
(i?s). The very broad line does not uniquely imply the 
presence of an accretion disk, however. Since the line 
does not possess a disklike, strongly double-peaked profile, 
other models may be able to explain the Ha profile. Spher- 
ical broad-line region (BLR) models incorporating radial 
outfiows can produce line profiles with horns and shoulders 
as well as very broad wings (Robinson, Perez, & Binette 
1990), qualitatively resembling the Ha profile in NGC 
4579. In addition, the broad-line profiles produced by the 
outflow models possess an intermediate-width core simi- 
lar in appearance to the "normal" broad-line component 
observed in NGC 4579. If the accretion disk model does 
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apply to NGC 4579, then the intermediate- width compo- 
nent of Ha most Ukely originates in the "normal" BLR of 
NGC 4579 rather than in the accretion disk, as the disk 
alone will produce a double-horned rather than centrally 
peaked profile (e.g., Eraclcous et al. 1995). Lincwidths of 
~ 2000 km s~^ are typical for the broad Ha components 
seen in ground-based spectra of LINERs (Ho et al. 1997b). 

Alternatively, models based on binary black holes 
(Gaskell 1983) can also produce very broad, asymmetric 
emission-line profiles. Such models have been ruled out, at 
least for a few well-studied objects, because the velocities 
of the broad wings do not exhibit the pattern of long-term 
variability that would be expected from orbital motion of 
a close pair of black holes (Eracleous et al. 1997). 

NGC 4579 shares other properties in common with the 
double-peaked emitters, most notably the overall shape of 
its spectral energy distribution (SED). In particular, NGC 
4579 lacks a "big blue bump" component, and its radio- 
to-optical flux ratio qualifies it as radio-loud although its 
radio luminosity on an absolute scale is small (Ho 1999b). 
Advection-dominated accretion flows (ADAFs; Narayan & 
Yi 1994) have been proposed as an explanation for these 
properties, and Quataert et al. (1999) have shown that 
the SED of NGC 4579 can be explained by emission from 
a thin accretion disk with a transition at ~ 100i?s to an 
ADAF for the innermost accretion flow. By truncating the 
thin-disk portion of the accretion flow and removing the 
innermost, hottest portion of the thin disk, the models are 
able to produce an SED without a big blue bump. Interior 
to the truncation radius, the disk is presumed to evapo- 
rate into an ADAF, and this portion of the accretion flow 
generates the X-ray emission. The broad Ha linewidth is 
consistent with this scenario, as it requires the truncation 
radius of the thin disk to be at r < 160 i?s. 

However, the value of Mbh used by Quataert et al. 
(1999) in modeling the SED of NGC 4579 is probably 
too low by an order of magnitude. The black-hole mass 
can be estimated using the recently discovered correlation 
between Mbh and bulge velocity dispersion (Ferrarese & 
Merritt 2000; Gebhardt et al. 2000). Central velocity dis- 
persion measurements for NGC 4579 range from 130 to 
187 km s~i (Whitmore, McElroy, & Tonry 1985; Gonzalez- 
Delgado & Perez 1996; Palacios et al. 1997; Heraudeau & 
Simien 1998). Applying the Gebhardt et al. (2000) and 
Ferrarese & Merritt (2000) relations, the predicted black- 
hole mass is in the range (1.8 — 9.3) x 10^ Mq.^ The smaller 
value of 4 X 10^ Mq used by Quataert et al. (1999) was 
a crude estimate based on the UV broad-line widths and 
ionization parameter (Barth et al. 1996), and the Mbh — cr 
correlations probably predict the central mass more accu- 
rately. With this larger central mass, the model predic- 
tions for the SED are altered, and it may be possible to 
fit the UV/optical portion of the SED with a thin disk 
extending all the way in to = 3i?s, as described by 
Quataert et al. (1999) for the case of M81. It remains to 
be seen whether models of a truncated thin disk surround- 
ing an ADAF can be constructed to simultaneously fit the 
UV/optical continuum and the X-ray luminosity of NGC 
4579, for a central mass of Mbh ~ 5 x lO'' Mq. 

Other evidence for an accretion disk in NGC 4579 is am- 



biguous. Iron K-line emission observed in ASCA spectra 
(Terashima et al. 1998, 2000) may originate from the inner 
accretion disk. However, the line does not clearly show a 
disk-like, relativistically broadened profile like those seen 
in some Scyfert 1 nuclei (e.g., Tanaka et al. 1995; Nandra 
et al. 1999), which would be expected if the disk extended 
all the way in to the black hole. Future X-ray observations 
with greater sensitivity and higher spectral resolution can 
better constrain the nature of the inner accretion flow and 
the possible presence of an ADAF. 

4. NARROW-LINE PROPERTIES 

Shock heating was the first mechanism proposed for the 
excitation of LINER narrow-line regions (NLRs), largely 
because observations of the [O III] AA4959, 5007 and A4363 
lines seemed to indicate electron temperatures higher than 
would be expected for photoionizcd gas (e.g., Koski & Os- 
terbrock 1976; Fosbury et al. 1978). Subsequent studies 
demonstrated that the observed [O III] line ratios could 
be reproduced by photoionization models provided that a 
high-density component (rie ^ 10^'^ cm~'^) is present in 
the NLR to boost emission in the [O III] A4363 line (Fil- 
ippenko & Halpern 1984). Correlations between narrow- 
line width and critical density for coUisional deexcitation 
("■crit), found in several LINERs, support this hypothesis 
by demonstrating the probable presence of radial density 
gradients in the NLR (Filippenko & Halpern 1984; Fil- 
ippenko 1985; Filippenko & Sargent 1988). Due to the 
difficulty of obtaining spectra with subarcsecond spatial 
resolution, however, there has been little direct evidence 
(in the form of spatially resolved density measurements) 
for density stratification in LINERs. The STIS observa- 
tions present an opportunity to revisit this issue, using the 
ratio of the [S II] AA6716, 6731 lines as a density diagnos- 
tic. 

Since the S/N of the observations is not sufficient for 
fitting the [S II] lines on a row- by-row basis, we performed 
multi-row extractions at several locations along the slit, 
for all five slit positions. When possible, we extracted 
regions corresponding to a well-defined cloud or clump. 
For each extraction region, we subtracted the local back- 
ground level and fitted the [S II] doublet using a pair of 
Gaussians. The two components were assumed to have 
the same velocity width, except at the nucleus where the 
higher S/N permitted a fit with fewer constraints. The 
line ratios were then converted to an electron density at 
each position using the STSDAS task TEMDEN, for an 
assumed temperature of = 10"* K. The [O I] A6300 fine 
intensity was also measured by fitting a single Gaussian. 
Due to the large difference between the critical densities of 
the [S II] (n„it = 1-6 X 10^ and 1.5 x 10* cm'^) and [O I] 
(j^crit = 1-6 X 10^ cm~^) transitions, the intensity ratio of 
[O l]/[S II] is sensitive to the nebular density as well. 

The results of these measurements are shown in Figure 
3. The electron density measured from the [S II] lines rises 
from < 200 cm~^ in the outer parts of the NLR (r > 100 
pc) to 2630 ±70 cm"-^ at the nucleus. Overall, the trend of 
increasing rie toward the nucleus directly supports the pic- 
ture of a radially stratified density in the NLR. Although 
the [S II] line ratio cannot probe densities greater than 
~ 10^ cm~^, the intensity and width of [O I] at the nu- 

®In applying the Gebhardt et al. (2000) relation, we neglect the small (< 10%) difference between the central velocity dispersion and the 
luminosity-weighted velocity dispersion within the effective radius, as spatially resolved measurements of a are not available. 
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cleus implies still greater densities in the inner NLR. The 
[O l]/[S II] ratio exhibits a sharp "spike" at the nucleus. 
Also, in the extraction centered on the nucleus, [O I] A6300 
(FWHM = 930 km s~^) is nearly twice as broad as the 
[S II] lines (FWHM = 500 and 570 km s^i for A6716 and 
A6731, respectively), indicating that the linewidth-ricrit re- 
lation continues to hold at the smallest radii resolvable by 
HST. The difference between the widths of the two [S II] 
lines in the nuclear extraction (an effect first seen in M81; 
Filippcnko & Sargent 1988) also demonstrates the pres- 
ence of a linewidth-ricrit correlation within the inner r = 2 
pc. 

A simple model for the density structure can be con- 
structed under the assumption that the ionization param- 
eter U (defined as the ratio of ionizing photon density 
to electron density at the illuminated face of a cloud) 
is roughly independent of radius (Filippenko & Halpern 
1984). Then, n,. oc r~^. The assumption of constant U 
is certainly an oversimplification, as the UV spectrum of 
NGC 4579 demonstrates that the NLR contains clouds 
with a range of values of U and rie (Barth et al. 1996), 
and the observed values of show considerable scatter. 
However, the model curve shown in Figure 3 demonstrates 
that an density law approximates the general trend of 
density as a function of radius. The [S III density in the 
nuclear extraction falls well below the r model predic- 
tion, but as noted already, the [S II] ratio at the nucleus 
most likely underestimates the central density. 

To confirm that the [O I]/[S II] ratio is indeed indica- 
tive of higher densities at the nucleus, we calculated a 
grid of photoionization models using Cloudy (Ferland et 
al. 1998) for a single, plane-parallel slab of solar-metallicity 
gas. Conditions were assiimcd to be typical of LINERs 
and of NGC 4579 in particular: an ionization parameter 
oi U = 10~^ — 10~^'^ and an ionizing continuum slope of 
aox = 0.9 (Barth et al. 1996). We find that [O I]/[S II] 
> 1 is only achieved at densities rig > 2 x 10** cm~'^ for 
dust-free gas, or at Ue > 10^ cm^^ when gas-phase de- 
pletion and dust grains with a Galactic dust-to-gas mass 
ratio are included in the calculations. Thus, the observed 
[O I]/[S II] intensity ratio of 1.01 ± 0.05 at the nucleus 
indicates that Ue continues to rise well above 10^ cm~^ 
within the unresolved, inner portion of the; NLR. and the 
inferred density increases further for the high- velocity ma- 
terial which contributes to the wings of the [O I] profile. 



5. CONCLUSIONS 

STIS observations have recently revealed very broad 
double-peaked or double-shouldered broad Ha profiles in 
three LINERs: NGC 4203 (Shields et al. 2000) and NGC 
4450 (Ho et al. 2000) from HST program GO-7361 (which 
observed a total of 24 targets), and NGC 4579 from pro- 
gram GO-7403 (with 6 targets observed). In contrast, no 
double-peaked emitters were found during the course of 
the Ho et al. (1997a) Palomar Observatory survey, which 
observed 486 nearby galaxies. Ground-based spectroscopy 
is apparently a very inefficient way to detect such features, 
except in much more luminous AGNs such as Arp 102B 
(Stauffer, Schild, & Keel 1983) and Pictor A (Halpern 
k Eracleous 1994). Among nearby {D < 20 Mpc), low- 
luminosity (Lboi < 10^3 erg s-i) AGNs, NGC 1097 is the 
only example of a double-peaked emitter discovered from 
the ground (Storchi-Bergmann et al.l993). It is difficult to 
determine meaningful statistics on the incidence of broad 
double-peaked emission in the LINER population, but an 
HST spectroscopic survey of a larger sample of galaxies 
would provide useful constraints. 

From the width of the very broad Ha line in NGC 4579, 
we infer an inner radius of ~ 160 Rg for the line-emitting 
portion of the accretion disk. Alternatively, BLR mod- 
els based on radial outflow rather than emission from a 
rotating disk may also be able to explain the observed 
Ha profile. Periodic HST monitoring of this galaxy (and 
other LINERs having double-peaked emission) would be 
valuable, as line-profile variability can provide tests of the 
accretion-disk and outflow models. 

Finally, the direct detection of a radial density gradient 
within the inner NLR of NGC 4579 provides support for 
the interpretation of LINERs as photoionized AGNs. 

Support for this work was provided through grant GO- 
7403 from STScI, which is operated by AURA, Inc., under 
NASA contract NAS 5-26555. NASA grant NAG 3-3556 
is also acknowledged. Research by A.J.B. is supported by 
a postdoctoral fellowship from the Harvard-Smithsonian 
Center for Astrophysics. We thank Gary Ferland for mak- 
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Fig. 1. — Spectra of the nucleus of NGC 4579. (a) STIS G750M spectrum through a 0'.'2 x 0'.'2 aperture. The Gaussian fits to the narrow 
lines and the "normal" broad Ha component arc shown. (6) The Ha profile after subtraction of the Gaussian narrow and broad components, 
(c) Ground-based spectra from Ho et al. (1997a), before and after starlight subtraction. 
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Fig. 2. — Comparison of the broad-line profiles of Ha, C IV A1549, and Ly a. The blue side of the Ly a profile is affected by Galactic 
absorption, and a spike of geocoronal emission has been removed from the spectrum. The narrow components of Ly a and C IV have not 
been removed. 
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Fig. 3. — Density diagnostics in the NLR of NGC 4579. (a) Electron density, determined from the [S II] line ratio, as a function of position 
along the slit. The solid curve represents the model rie oc r~^. (b) Intensity ratio of [O I] A6.300 to [S II] AA6716, 6731. In both panels, the 
horizontal error bars represent the extraction width used for the emission-line measurement. Vertical error bars represent the uncertainties 
in the measured emission-line intensities. At the assumed distance of 21 Mpc, 1" corresponds to 101 pc. 



